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Abstract 

The evolution of QCD coupling constant at finite temperature is con- 
sidered by making use of the finite temperature renormalization group 
equation up to the one-loop order in the background field method with 
the Feynman gauge and the imaginary time formalism. The results are 
compared with the ones obtained in the literature. We point out, in par- 
ticular, the origin of the discrepancies between different calculations, such 
as the choice of gauge, the break-down of Lorentz invariance, imaginary 
versus real time formalism and the applicability of the Ward identities at 
finite temperature. 



1 Introduction 



One of the application of field theory at finite temperature [1-4] is to find the behaviour 
of coupling constant as a function of energy, temperature, and the chemical potential 
using the renormalization group (RG) equation. The knowledge of coupling constant 
at finite temperature environment then can be used, for instance, in perturbative 
calculations for quark-gluon plasma created in ion-ion collisions at high energies, in 
the evaluation of the grand unification scale in a cosmological context and the shift 
of the energy levels in a hydrogen plasma. 

Recently a great deal of attention has been paid to the study of the behaviour 
of the coupling contant at finite T [5-19]. The resulting formula for the temperature 
and scale dependent part of the coupling constant has proved to be very sensitive to 
the prescription chosen. In this paper we wish to examine this and clarify the origin 
of the discrepancies between different calculations in the literature. We employ the 
background field method [20,21], which is based on a manifestly gauge invariant gen- 
erating functional, and the Feynman gauge. The background field method provides 
notorious simplifications since we have to calculate only the renormalization constant 
of the background field gluon propagator. We discuss this method in section 2.1. 
Since the coupling constant depends not only on energy but also on temperature, we 
derive a pair of RG equations - one for energy and one for temperature. This is done 
in section 2.2 and in section 2.3 we derive its solution. In section 2.4 we discuss the 
Feynman rules at finite temperature in the imaginary time formalism. 

In order to get the renormalization constant of the background field gluon prop- 
agator we have to calculate the polarization tensor. At finite temperature due to the 
lack of the Lorentz invariance the structure of the polarization tensor is not equiva- 
lent to the zero temperature one. In section 3.1 we give a prescription how to define 
the renormalization constant at finite temperature. In section 3.2 we calculate the 
polarization tensor in the one-loop approximation. In section 3.3 calculating the vac- 
uum part of the polarization tensor we reproduce the well known zero temperature 
formula, i.e. the standard vacuum QCD one-loop result [22,23], which reads as 

[9\t^)V = k'(/^o)]-Ml + 25o'(/^o)i^o/(47r)2] , (1.1) 

where 

Ko = (117V/3 - 2nf /3) ln(/////o). (1.2) 

Here N is for SU{N), is the number of flavours, ji is the energy scale and jiQ is 
the reference scale. 

We calculate the matter part of the polarization tensor in section 4.1 and, in 
section 4.2, we derive the temperature and scale dependent part of the coupling 
constant. 

In section 5.1 we review the results of the previous works and in 5.2 we compare 
the asymptotic expansion formulas &X, a — ^/T <^ 1 derived in different schemes. 
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We give a few comments and, in particular, discuss the origin of the discrepancies 
between different calculations. 

In Appendix we collect the relevant integrals which we encounter in the text. 



2 Formal methods 

2.1 Background field method 

Explicit gauge invariance, which is present at the classical level in gauge field theories, 
is normally lost at the quantum level. This can be seen from the generating functional 
for QCD: 

Z(J) = J DAD^D^Dr]Dfiex.p{i J d'^x[jCQCD(A,^f, 

-(G'^)V(20 + J^'^A.a + ndGydoj'^r)]}, (2.1) 

where A is the gluon field, ^' and ^ are the fermion fields, rj and f] are ghost fields, 
(G")^ is a gauge fixing term (^:gauge parameter), a;" is a SU{3) group parameter, 
and J'^A^a is a source term. In SU (3) the gauge fields transform in the following 
way: 

A^^U{A^-idJg)U\ (2.2) 

^U-^,^ (2.3) 

r]^Urj,fj ^fjU\ (2.4) 

where U{x) = exp[icj(a;)"Ta] is a unitary transformation, and Ta is a generator for 
SU{3). By using transformations (2.2)-(2.4) one can see that the gauge invariance of 
Z{J) is lost in commonly used gauges such as d^A^j^ = 0. 

The advantage of the background field method [20,21] is that it can maintain the 
explicit gauge invariance. For this purpose we devide the gluon field ^4^ into a sum 
of a classical background field and a quantum field Q/^ 

A^^B^ + Q^, (2.5) 

and choose for the gauge fixing term G" the background field gauge condition 

G« = {D,{B)r''Ql (2.6) 

where is the covariant derivative: 

{D^,)ab^d^5,f, + gUcB;. (2.7) 

The generating functional reads 

Z{J,B) = I DQD^D^DriDf] J exp{i J d^'xlCgcDiB + Q,^ 
-{G^f/m + J:AI + ndC^/du^rj]}. (2.8) 
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This functional is gauge invariant, which follows from the transformations: 



B^^U{B^-idJg)U\ (2.9) 

^ UQ^U\ (2.10) 

D^{B) ^ UD^{B)U\ (2.11) 

J^{B) ^ UJ,{B)U\ (2.12) 



together with Eqs. (2.3) and (2.4). Notice that we are only changing the integration 
variables in Eq. (2.8). 

In the background field method the quantum gauge fields and the ghost field need 
not be renormalizcd since they appear only inside loops. No vertex functions are to 
be considered. Thus only renormalizations 

9o = ZgQR, (2.13) 
Bo = {ZbY'^Br, (2.14) 
6 = Z^U, (2.15) 

are needed. 

The explicit gauge invariance of Z{J,B) implies that the perturbation series is 
gauge invariant in every order in and that the background field renormalization 
factor Zb and the coupling constant renormalization Zg are related with each other. 
The field tensor F^^, in jCqcd has to be gauge covariant and is renormalized as 

= (ZsY/'mBtU - d.iB^U + 9RZ,iZB f'r'^iB^)RiB;)n]. (2.16) 
Thus we have a relation 

Zg = {Zb)-^I\ (2.17) 

which enables us to calculate the evolution of coupling constant g^i in a simple way. 
We assume this relation to be valid also at finite temperature. 



2.2 Derivation of the RG equations 

Prom Eqs. (2.13) and (2.17) we have 

% = 9r{.Zb)-^I''. (2.18) 

We use the dimensional regularization [24,23] and perform our calculations in (4 — 2s) 
dimensions. We notice from the Lagrangian of QCD that the dimension of a gluon 
(quark) field is , where /i is the scale parameter. We redefine go so that 

it becomes dimensionless and rewrite Eq. (2.18) as 

9o^9r{Zb)-'^^I^', (2.19) 
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where the bare couphng qq does not depend on temperature T and the scale /j,. A 
pair of RG equations results upon taking the derivative of Eq. (2.19) with respect to 
T and fj, 

Td/dT[gRZs'^^fj,'] = 0, (2.20) 
fxd/dfxigRZ-'^^,,'] = 0. (2.21) 

These equations, which are not symmetrical in T and /i, determine the behaviour of 
coupling constant qr with respect to T and /i changes. Generally Zb has the form 
[22,23] 

oo i 

= 1 + E E + f^Ht^, n, (2.22) 

i=l j=l 

where je^ are the divergent contributions independent on T and while /'^•'-'(/U., T) 
are convergent temperature and scale dependent functions which vanish in the low- 
temperaturc limit. 

From Eqs. (2.21) and (2.22) one can derive the following formula, by taking e — > 0, 
(from now on we write g — gR, for simplicity) 

l^dg/di, = g\-A^'^ + (/./2)a/a/.[/(^)(/., T)]} + 0(/). (2.23) 

This equation is reduced to the ordinary one-loop RG equation in the zero- temper- 
ature limit since in this limit the function /*^^^(/i, T = 0) vanishes. 

Similarly one can derive the RG equation for T from Eqs. (2.20) and (2.22) 

Tdg/dT ^{g^/2)T^/^T[f^^\^Ji,T)\ + 0{g'>). (2.24) 

Equations (2.23) and (2.24) constitute the coupled RG equations. 



2.3 Solution of the coupled RG equations 

The solution to Eqs. (2.23) and (2.24) can be found in a straightforward way by 
integrating Eq. (2.23) from jiQ to ji and Eq. (2.24) from Tq to T. We have 

l/[2/(/xo,T)] - l/[2/(^,T)] = ln(/.//xo) 

+ [/(^)(/.,T)-/W(/xo,T)]/2, (2.25) 

l/[2/(/x,To)] - l/[2/(^,T)] = [f^'\^^,T) - /«(/i, To)]/2. (2.26) 

Putting T equal to Tq in Eq. (2.25), and using Eq. (2.26) we arrive at the desired 
solution 

[g\li,T)]-^ = [gHf^o,To)]-' + 2A(i) ln(/.//.o) 

-[f^'K^^,T)-fi'\^^o,To)]. (2.27) 
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This equation describes the evolution of the QCD couphng constant as a function 
of the momentum scale and the arbitrary temperature [5-9]. (/io,2o) denotes the 
reference point. 



2.4 Imaginary time formalism 

The Feynman rules at finite temperature T{— are derived from the generating 
functional 

Z{J) ^cj DQD^D^DrjDfiex.p[J^ dr J d^x{C{x,T) 

-[G''ix,T)]yi20+vix,T)dGydLu''vix,T)}], (2.28) 

where r], fj, Q are periodic fields, while \1', ^ are antiperiodic fields. 

This functional differs from Eq. (2.8) only in that time r is now imaginary it. 
Infinite time domain has been compactified to the finite interval [0,(3]. Prom Eq. 

(2.28) one can formulate the Feynman rules at finite temperature by modifying the 
ones at zero temperature in the following way [25] : For the loop integral we have 

/ d'p/(2n)' ^ (i/P) E / dV(27r)^ (2.29) 

and for the time-component of 4-momentum we have for bosons and ghosts: 

Po = i2n7r//?, (2.30) 

and for fermions: 

Po = «(2n + l)7r//3 + ii^n, (2.31) 

where ^ch is a chemical potential. The frequency sums for bosons and fermions in Eq. 

(2.29) are readily changed to contour integrals [26]: We have for bosons and ghosts 

{i/P) ^ /(po = ^2r^7r//3) = (l/27r) / dpof{po)NB{po) 

n=-oo J-ioo+e 

/200 — E pioC 
dpof{po)NB{-po) + (l/27r) / dpof{po), (2.32) 
-ioo—e J—ioo 



where 

and for fermions 



7VB(po) = l/[exp(/?po)-l], (2.33) 

oo 

(i/P) E /[Po = ^(2r^ + l)7^//3 + /.eft] = -(l/27^) 

n=— oo 

/ioo+e rioo—e 
dpofipo + f^ch)NF{po) + (l/27r) / dpo 
-joo+e J —ioo—£ 

/ioo 
dpofiPo), (2.34) 
-ioo 
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where 

iV^(po) = l/[exp(/3po) + l]. (2.35) 

The first two terms in the right-hand side of Eq. (2.34) correspond to particle and 
anti-particle contributions respectively and vanish at T = 0. Since we will be only 
interested in phenomena in an enviroment where the sum of quantum numbers is not 
conserved, hereafter we put fich equals to zero. Also we have omitted the term which 
gives the finite density contribution at T = 0. 



3 Polarization tensor 

3.1 Structure of the polarization tensor 

At a zero temperature environment the polarization tensor is Lorentz invariant and 
can be expressed as [22] 

n«t = 5"'n^., (3.1) 

where 

n^, = n(A;V-W- (3-2) 

The zero temperature polarization tensor is transverse with respect to k (current 
conservation) : 

A;^n^. = 0. (3.3) 

At finite temperature, in the presence of matter, the Lorentz invariance is lost and 
the polarization tensor can only be 0(3) rotational invariant [27,17]. Then the polar- 
ization tensor can generally depend only on 4 independent quantities, which we can 
choose, for example, IIoo, kjloi and the two scalars 11^ and 11^ appearing in 

Uij = nr(% - kikj/k'') + ULkikj/k^. (3.4) 

At finite temperature whether the transversality condition is satisfied or not depends 
on the gauge used. The polarization tensor is not transversal, e.g., in the Coulomb 
gauge {diAi = 0), but is transversal at the one- loop level in the temporal axial gauge 
{Aq = 0), and in every order of the perturbative calculations in the background field 
gauge [17]. The transversahty condition (3.3) restricts the structure of the polariza- 
tion tensor. Prom Eq. (3.3) we have 

koUoo = kiUio, (3.5) 

and 

klUoo = kikjUij. (3.6) 
Using Eqs. (3.4) and (3.6) we obtain 

Hl = klUoo/k\ (3.7) 
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For the coefficient of the transversal part of the polarization tensor we get from Eqs. 
(3.4) and (3.7) an expression 

Ut = (Uu - klUoo/\^^)/2. (3.8) 

By using 11^ = (H^^ — IIoo), we rewrite Eq. (3.8) as 

Ut = [n^^ - noo(l + kl/k^)]/2. (3.9) 

Thus in Eqs. (3.4)-(3.9) we have derived the general form of the temperature depen- 
dent 0(3) symmetric polarization tensor in the background field method. 

The polarization tensor can be splitted into a sum of a temperature independent 
(vacuum) part and a temperature and scale dependent (matter) part [4] : 

n^;/|^(jc ~i~ n^;/|jjja^(. (3.10) 

At zero temperature limit the temperature dependent matter part vanishes. 

The structure of the polarization tensor at zero temperature was given in Eqs. 
(3.1) and (3.2). The polarization tensor is related to the gluon propagator, D^,^, as 

= - Doi. = ^D^U- (3.11) 

The relation between the bare and the renormalized propagators (see Eq. (2.14)) is 

Dqhv — Z^^Dr/j,!,, (3-12) 

which together with Eq. (3.11) leads us to 

Zs = 1 - n. (3.13) 

Thus the renormalization constant Zb can be obtained from the (background field) 
gluon self-energy tensor. 

In investigating the behaviour of the coupling constant one has to include the tem- 
perature dependent parts of the polarization tensor in the renormalization constant 
Zb via Eq. (3.13). Although the finite renormalization is somehow arbitrary one has 
to follow some rule or prescription in every order of perturbation to be consistent 



At finite temperature we encounter another type of ambiguity which is caused by 
the lack of the Lorentz invariance. In order to define the polarization tensor at finite 
temperature we generalize Eq.(3.13) for T = as: 

Zb — i — ^\vac — ^\matt, (3-14) 

where we have either 

nUatt = Hoo/k^, (3.15) 

or 

nUa« = nr/kl (3.16) 

Naturally at T = 0, Eq.(3.14) reduces to Eq.(3.13). IIoo and IIt are not connected 
with each other and in general there is no a priori way to decide which one is more 
natural. 
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3.2 The polarization tensor at the one-loop level 

Our working formulas, which allow to analyze the evolution of QCD coupling constant 
at finite temperature, are:Eq.(2.22), Eq.(2.27) and Eqs. (3.14)-(3.16) for T ^ case 
[Eq.(3.13) for T = case]. Accordingly we have to evaluate the self-energy diagrams 
in Fig.l to get the renormalization factor for Zb up to the one- loop order. 

The Feynman rules for the interaction vertices are the same as in the zero tem- 
perature case, and therefore identical to those given in [21]. Evaluating the one-loop 
diagrams (la)-(ld) in the Feynman gauge (i.e., we set ^ = 1 in the gluon propagator) 
we find for the boson contributions in the polarization tensor 

^^v\hosan = ig^^ j ciV(27r)^[%^A;^ + 2{k^p^ + Kp^) 

+Ap^p, - Zk^K - 2{k + pfg^MD, (3.17) 

where 

D={k + pfp^. (3.18) 

For the fermion loop, neglecting the quark masses compared to momentum scale and 
temperature, we have from the diagram (le): 

^Hu\ fermion = -^ig^TpUF j d^p/ {2n)^{k^p,, 

+k,p^ + 2p^p, - g^,[{kp) + p'W/D, (3.19) 

where Tp = 1/2 for SU{3). We have also included np (number of fiavours) in Eq. 
(3.19). 

Our polarization tensor satisfies the transversality condition 

^ljJ^lxv\boson k/jJ^iMul fermion 0- (3.20) 

This can be shown exphcitly as follows: From Eq. (3.17) we have 

Til^Aboson^ig'N J d*p/{2TTY[-kJp' 

-ik+ p)J{k + p)2 + 2p,/p2 - pj{k +pf]. (3.21) 
Replacing k + p ^ p in the second term we have 

^^.^boson^-^g''N j d''p/{2T:f[{k-p)Jp'+pJ{k+pf]. (3.22) 

Then by changing k — p ^ —p in the first term, we immediately come to Eq. (3.20) 
for the boson contributions. Similarly one can prove the transversahty condition for 
the fermion polarization tensor. In the next subsection we will extract from 11^,^ the 
information on the vacuum part of the polarization tensor. 
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3.3 Vacuum part of the polarization tensor 

Using a standard technique and introducing the Feynman parametrization in Eqs. 
(3.17) and (3.19), namely 

1 /(ah) = C dx{l/\ax + 6(1 - x)f}, (3.23) 

we have 

l/D = l/\p^{k + pf] = C dx{l/\(p + kxf + k'^xil - x)?}. (3.24) 

JO 

Changing variable p ^ p—kx we get the following formula for the boson contributions 
from Eq. (3.17): 

^nulboson = ig'^N J dV(27r)^ dx[Ag^,e + 2(1 - 2x){kf,p, + Kp^) 
+^P^.Pu + (-3 -4x + ix^)k^k,]/{p^ + K^f, (3.25) 

where 

= k'^x{l- x). (3.26) 
Similarly we get the formula for the fermion contributions from Eq. (3.19): 

^nA fermion = -^ig'^up J (fp/{27rY dx[2{-x + x^)k^ki, + (1 - 2x) 

y<{k^Pu + Kp^) - g^v{kp) + 2p^p, + g^.k^xl/ip" + K^f. (3.27) 

In these equations for the vacuum parts wc have dropped the terms proportional to 
1/p^ and 1/ {k + pY, which turn out to be zero after simple calculation [see Eq. (A.l) 
in the Appendix]. Notice that all of the integrals are ultraviolet divergent and thus 
have to be regularized. For this purpose we employ the dimensional regularization 
[24,23], which preserves gauge symmetries explicitly. The integrals for the vacuum 
parts become Euclidean if we change ipo Pa and thus can be easily evaluated. We 
obtain the following results: 

k^k,)/e + 0{l), (3.28) 

k^k^)/e + 0{l), (3.29) 

(3.30) 

Combining Eqs. (3.28) and (3.29) and taking into account Eqs.(2.22), (3.2) and (3.13), 
we get for A^^) in Eq. (2.27) as: 

>lW ^ (117V - 2nF)/(3Q;), (3.31) 



^ixAboson = -'^'^g'^N/{2,a){g^^k'^ - 

and 

fermion = 2g'^nF / {'ia){g^„k'^ - 

where 

a = (47r)2. 
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which is, as expected, in accord with Eq. (1.2) for zero temperature. In concluding 
this subsection let us emphasize the simplicity of the calculation by the background 
field method in contrast to the conventional methods, e.g., in the covariant gauge 
[7-9]. 



4 Temperature dependent parts of the polariza- 
tion tensor 

4.1 Calculation of IIoo and 11//// 

In this subsection wc calculate IIoo and H^^, without using the Feynman parametriza- 
tion. To specify the subtraction point we employ the static limit of zero external en- 
ergy, which is commonly used in the literature [5] . In this prescription the momentum 
k is specified to be space-like k = {0,ki){i = 1,2,3) with k"^ = — /x^. Such a choice 
enables us to determine the static properties. Using Eqs. (3.17) and (3.19) and the 
integrals from the Appendix we obtain for IIoo 

Hoo = 2tg'N J d'p/{27ry[2k' + 2pl - (p + ky]/D\,,,^ 

+2ig^nF J d'^p/{2T:Y[-2pl + {kp)+p^]/D\ 

fermion 

= g^T\N/6 + np/U) - 2g'' N ii''[AFBo(a) - ^^2(0)] 
-2g^nFAFFo{a) - Fp^ia)], (4.1) 

with a — (3/1. Similarly we derive an expression for 11^^: 

= W'N j d^p/{27:f[^e - 12{kp) - V]/L>|,„^„„ 

^Aig^np j d'^p/ {2Tif \{kp)+p']/D\ 

fermion 

= g''T\N/'i + np/Q) - 2^V[ll^^^so(a) + 2n^F^o(a)], (4.2) 
where we define the boson functions 

roo 

FsniPpi) = (l/a) / dxx''NB{pLx/2)L, (4.3) 
Jo 

with 

L = ln|(l + x)/(l-x)|. (4.4) 

The integrals (4.3) are not ultraviolet divergent and hence do not give infinite 1/e 
contribution. The fermion functions Fpnio) are defined by replacing in Eq. (4.3) 

iVs(W2) ^ NF{f,x/2). (4.5) 
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4.2 Temperature dependent part of the coupling constant 



As pointed out in section 3.1 we encounter an ambiguity in determining the renormal- 
ization constant as a direct consequence of the lack of Lorentz invariance. Here 
we write two formulas, one derived from IIoo and the other one from n^-. Prom Eq. 
(3.15) we have (a = 

-f^^\\i, T) = {N/6 + nF/12)/a^ - 2iV[4Fso(a) - Fs2(a)] 
-2nF[FFo{a) - FF2{a)], (4.6) 

and the other one from Eq.(3.16) and 11^ = {^nn — noo)/2 

-f^'\l2, T) = (iV/12 + nF/2A)/a^ - N[IFBo{a) + ^52(0)] 
-nF[FF^{a) + FF2{a)]. (4.7) 

The result for Eq. (4.7) was reported in a short communication in [18] (where an 
overall factor of 1/2 was missing in the right-hand side of formulas (8) and (9)). 



5 Discussion 

5.1 Review on the results of previous works 

Gendenshtein [5] calculated the QCD coupling constant at finite temperature in the 
one-loop approximation by using the RG equation, the dimensional regularization, 
and the covariant gauge with a space-like normalization momentum — (0, //) and 
obtained 

A(^) = (llAr-2nF)/(3a), (5.1) 

as in Eq. (3.31), and 

-f^\^^,T) = {N/^ + nF/Q)/a\ (5.2) 

where we have presented the results using our notations. 

Kajantie et al. [6] studied the gauge field part of QCD with N colors. They used 
the Aq — Q gauge and defined two renormalization schemes by writing 

D-l = D^^^{Za - G/p') - (F - G)Pi., (5.3) 

where D^i, is the gluon propagator and F, G and come from the polarization 
tensor 

n^. = FPj:, + GPl. (5.4) 

In the "magnetic prescription" they fixed the propagator at the point p^ = (0, n) and 
had for the temperature depending function the expansions (without quarks) 

-/W(//,r) = 5A^/(16a), fora<l, (5.5) 
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-/(^)(/x,T) = Ar[l7/(90a^) + 83a/(6300a^)], for a > 1. (5.6) 
In the so-called "electric prescription" they derived the couphng constant from F as 

-/(^)(/x,T) = iV[l/(3a^) - l/(4a) -22(lna)/(3a)], for a < 1 (5.7) 

- f^\i^,T) = -N[l/{18a^) - lla/(900a^)], for a > 1. (5.8) 



Notice here the change of sign in the high T behaviour. They concluded that the 

"magnetic prescription" is more natural than the "electric" one because the former 
uses the physical part of the gluon propagator. They also argued that one has to 
choose iJ, = 3T in the thermal equilibrium. In that case the effect of the temperature 
dependent parts becomes negligible. 

Nakaggawa, Niegawa and Yokota [7] used the real-time formalism and studied the 
scale-parameter ratios A(a)/A derived from different vertices in 4-flavour QCD. They 
used the covariant gauge and found that the ratios derived from the three-gluon vertex 
and the gluon-ghost vertex show just the opposite behaviour than the one derived 
from the gluon-quark vertex. Namely in the former case one gets a growing ratio 
while in the latter case a decreasing ratio. 

Fujimoto and Yamada [8] used the real-time formalism and derived the tempera- 
ture depending coupling constant from the gauge-independent Wilson loop. At a -C 1 
, it reads as (without quark contributions) 

- /«(/., T) = C[l/(3a^) - l/(4a) - 22(lna)/(3a)]. (5.9) 

The same authors have discussed the finite temperature RG equations [9] in the 
one-loop approximation, using the real-time formalism, the covariant gauge and the 
dimensional regularization. Prom the gluon propagator and three-gluon vertex they 
obtained: 

-/(i)(/i, T) = (C + Tf)/{4a^) + C[-23/3FBo(a) - 3FB2{a) - U/3GBo{a) 
+32^^2(0)] + 2Tf[FFo{a) - SFp^ia) + 152/9G'fo(a)], (5.10) 

where Tf = 1/2, and C = for SU{N). The coefficient in front of FBo{a) was 
originally 4/3. The error was pointed out in Ref. [10]. 
Prom the fermion propagator and fermion-gluon vertex: 

-/«(/i, T) = {C + 3Cf + 2Tf)/{12a^) - 2{Cf + llC)FBo(a)/3 
-16{Cf + 10C)GBo(a)/9 - 32CGB2{a) + 2(-5C/6 + C//3 + Tf) 
xFpoia) - 4{SCf + 17C)GFo{a)/9 + 16CGF2{a), (5.11) 

where C/ = (A^^ _ i)/(2N) for SU(N). 



12 



From ghost propagator and ghost-gluon vertex: 



T) = (C + 7»/(12a2) - C[7FBo{a) + FB2{a) 
+2GBo{a)] - 2Tf[FFo{a) + FF2{a)]. 



(5.12) 



In the above equations the boson functions FBn{ci) are defined in Eq.(4.3), while 
GBn{o) are defined as 



Fcrmion functions Gpnio) are defined by replacing NB{^Jix/2) by Np{nx/2) in Eq. 
(5.13) as FBn{ci) was defined from Eq. (4.3). The functions Gin{A){i = B,F) do not 
appear in our results Eqs. (4.6) and (4.7), since they come from trigluon renormal- 
ization only. [Note an amusing coincidence: -J^^) in Eq. (4.6)= -/(^)[G'Bo(a) ^ 0] 
of Eq. (5.12), both of which are derived from very different prescriptions.] 

Stephens ct al. [19] performed a background field one-loop calculation of gauge 
invariant beta functions at finite temperature, using the retarded/advanced formalism 
developped by Aurenche and Becherrawy [28]. In terms of our notations, their result 
reads as: 



The numerical coefficients of the leading terms of the high-tempcraturc expansion 
and the fermion parts are in complete agreement with our result in Eq.(4.7). In the 
boson parts, however, there are some numerical discrepancies with our result. 



5.2 Comparison of asymptotic expansions 

In order to compare our results with those mentioned in section 5.1 we derive the 
asymptotic expansions for —f^^\fi,T). The asymptotic expansions for Eqs. (5.10)- 
(5.12) at a = /i/T <S 1 in the high-temperature regime read as follows: 



GBn{a) = {2/a) / dx / dyx^+^NB{i^x/2)/[x'{y' + 3) - 1] 
Jo Jo 



(5.13) 



-/(^) [fi, T) = {N/12 + nF/24)/a^ - iV[21/4Fso(a) 
+FB2{a) + 7/2GBi{a)] - n^[Fj.o(a) + ^^2(0)]- 



(5.14) 



-fW^^^T) = 7iC/{9V3a^) - 23C/(48a) + [(-22 
+7r/3V3)G + (18 - 387r/9V3)T/](lna)/(3a), 



(5.15) 



/«(/i,T) = [(l-7r/v^)C + 3C/]/(12a')-(llC + C/)/(24a) 
[(23 + 827r/9V3)C + (4 - 87r/9V3)C/ + 4r/](lna)/(3Q;), 



(5.16) 



/W(//,r) = -7C/(16a)-[C(22 + 7r, 



yVS) -87>](lna)/(3Q;). 



(5.17) 
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Next we derive the asymptotic expansions for our results, i.e., Eqs. (4.6) and (4.7) at 
a <C 1, by using the following high temperature expansions [9]: 

FBo(a) = l/(16a) + (£i - l)/a, (5.18) 

FbM = l/(12a2) + (£i - l/6)/(3a), (5.19) 

Fpoia) = -{i2 - l/2)/«, (5.20) 

FF2{a) = l/(24a2) - (£2 - l/6)/(3a), (5.21) 

where 

ii = ln(a/47r) + 7, (5.22) 

4 = ln(a/7r) + 7. (5.23) 
We find the asymptotic formula from Eq. (4.6) 

-/(^) (/i, T) = {N/3 + n^/6)/a^ - iV/(2a) 

-[N{22£i - 71/3) - 4n/(£2 - 2/3)]/(3a). (5.24) 

The term coincides with the one of Gendenshtein [5], i.e.,Eq. (5.1). The gluon 
part in this asymptotic formula is consistent with the result of Elze et al. [17] derived 
in the background field method, which in our notations reads 

-/W(//,r) = A^[l/(3a2)-l/(2a)-22(lna)/(3Q;) + ...], for a < 1. (5.25) 

This coincides with the result of Nadkharni [12], who has the terms up to 0(1 /a). 
The dehcate reason why in the background field method we have a factor -1/2 in 
front of the second term proportional to 1/a, while in some results [see Eqs. (5.7) 
and (5.9)] it is equal to -1/4 is clarified by Elze et al. in [17]. 
Next from Eq. (4.7) we find the asymptotic formula: 

- /W(//, T) = -7A^/(16a) - [N(22ei - 127/6) - 4/1^(^2 - 5/12)]/(3q;). (5.26) 

The behaviour of the transverse polarization tensor IIt in the infrared region is 
known to have a form [17] 

limnT(0,k)//x2 - g'^[-cN/a + 0{\na)]. (5.27) 

The factor c has been calculated in different gauges. In the covariant ^-gauge its value 
is [27] 

c=(9 + 2e + a/64, (5.28) 
whereas in the temporal axial gauge it is [6] 

c = 5/16. (5.29) 
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Our formula (5.26) at small a behaves as Eq. (5.27) with 



c — — 



7/16, for bosons, 



(5.30) 



and 



c — 0, for fermions. 



(5.31) 



These results can be read also from the results of Refs. [9] [see Eq. (5.17)] and [12]. 
The relation (5.31) was also noticed by Elze et al. in [17]. We remark here that our 
c is negative and hence no spurious pole appears in the transversal propagator 



in contrast to the covariant and the temporal axial gauge cases. 

It is a known fact that the terms are gauge independent and the gauge param- 
eter dependence starts at ~ T order [13]. Thus the fact that the terms in Eqs. 
(5.15) and (5.16) are different from others, e.g. Eq. (5.1) should not originate from 
the gauge choices. 



5.3 Concluding remarks 

In concluding the paper we give a few comments: 

1) In sections 5.1 and 5.2 we have seen that the temperature dependence of the 
QCD coupling constant is very sensitive to the prescription chosen. This is not a 
trivial issue, because all the results obtained hitherto also heavily depends on the 
vertex chosen (i.e., the trigluon, the ghost-gluon, or the quark-gluon vertex) to satisfy 
the renormalization condition of the QCD coupling constant. Furthermore in some 
gauges, e.g., in the Coulomb gauge, the transversality condition (3.3) on the polariza- 
tion tensor does not hold and hence the structure of the polarization tensor in such 
gauges is different from that which satisfies the condition. 

One of the reasons why we encounter different results in the literature is that 
the broken Lorcntz invariance has not been treated as wc have done in Eqs. (3.15) 
and (3.16). For example in deriving Eqs. (5.10)-(5.12) the authors of Ref.[9] have 
extracted the gluon and the vertex renormalization constants in front of Lorentz 
invariant structures not paying attention to the break down of Lorentz invariance. 

Another possible source for discrepancies is that the Ward identities are used 
at finite temperature for different renormalization constants. Remember that the 
derivation of the Ward identity is based on the gauge invariance and also on the 
Lorentz invariance at T = 0. Thus it is not a surprise that results from different 
gauges or even (within the same gauge) from different vertices are totally different. 

The correspondence between the imaginary and real time formalisms has been 
investigated in detail [14-16]. The choice of the formahsm can also be the source of 
the discrepancies under study was pointed out in these references. 



DT^,(o,k) = -{s,, - kkj/k^)/[k^ - nT(o,k)], 



(5.32) 
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To clarify the issue of choosing a suitable renormalization prescription, one would 
need to compute the two-loop contributions to the coupling constant at finite tem- 
perature. It was shown, in particular, in the massive 0{N) scalar model that the 
one-loop result is drastically changed by two-loop contributions at high T and in zero 
momentum hmit [29]. 

2) As we have seen the function f^^^ {/j,, T) shows a power-like T-dependence instead 
of a logarithmic fall-off as a function of /x. Thus we have 

= E c.[(r//.r-(ro///or] + .... (5.33) 

n=l,2 

However, if a relation T = {const) fi holds, where (const) ~ 1/3 [6], then one would 
have a logarithmic T-dependence like in a zero temperature environment. Such a 
relation holds if one is dealing with the thermal equilibrium. In a thermostat one 
cannot speak about individuals and thus n is not a measurable quantity. 

We have derived Eq. (2.27) using two essentially different RG equations (2.23) and 
(2.24). Suppose we had used only one RG equation for T, i.e., Eq. (2.24), then we 
would have instead of Eq. (5.33) 

-[f'\^l,T) - f'\^lo,T)] 

= E c„[(T//.)"-(r//.o)1 + .... (5.34) 

n=l,2 

In this case we would have a power-like dependence even with the relation T — 
(const) II. 

3) Since in a thermostat T is a measurable quantity, though fi is not as, stated 
in 2), we should be able to define the coupling constant as a function of only T. 
This could be done by calculating an expectation value of the coupling constant. In 
that case as a probability density function we could use either the Bose or the Fermi 
distribution. However such an expectation value of the coupling constant could not 
be used to a scattering process in a thermostat wherein a particle with some definite 
energy enters. 

We would like also to mention that all the machinary for the evolution of the 
running coupling constant at finite temperature can be analogously applied for the 
case of a quantum field thcory(such as QCD) at finite energy as formulated in [30]. 

In conclusion, we note in accordance with the previous observations that there 
is in fact no unique way to define a temperature depending QCD coupling constant 
and the issue of finding its reasonable prescription is left as a subject of further 
investigation. 
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Appendix 



The formulas for the integrals which we encounter in the text are assembled in 
this Appendix. 

Jd'p/{2nr[l/p']U = 0, (A.l) 

J d'p/{2nni/p%oson = ^TVl2, (A.2) 

I d'p/{27r)'[l/p%ermion = -iTy2A, (A.3) 

with Hch = 0. 

In the following we use the static space-like prescription specified by /cq = and 

I dV(27r)'(l/L')|6o.on = -2iFBo(a), (A.4) 
where FBo{a) is defined in Eq. (4.3). 

J d''p/{27r)^[{kp)/D]\,oson = -i//'FBo(a), (A.5) 

J d'p/{2Ti)\pl/D)\,,,^ = -z(/.V2)Fso(a). (A.6) 

For the fermions we have similar results as in Eqs. (A.4), (A.5) and (A.6) in which 
we just replace Nb{ij>x/2) by Nf{iix/2) and change the sign in front of the equations. 
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Figure Caption 

Fig.l. Diagrams for the one-loop calculation of the background field renormalization 
factor Z^'- a) Gluon loop; b) Ghost loop; c) Gluon loop from the 4-gluon vertex; 
d) Ghost loop from the 2-gluon 2-ghost vertex; e) Fermion loop. Wavy lines are 
quantum gauge field propagators. Wavy lines terminating in a " B" represent external 
gauge particles. Solid lines are fermion propagators and dashed lines represent ghost 
propagators. 
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